Accuracy vs run time in adiabatic quantum search 
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Adiabatic quantum algorithms are characterized by their run time and accuracy. The relation between the 
two is essential for quantifying adiabatic algorithmic performance, yet is often poorly understood. We study 
the dynamics of a continuous time, adiabatic quantum search algorithm, and find rigorous results relating the 
accuracy and the run time. Proceeding with estimates, we show that under fairly general circumstances the 
adiabatic algorithmic error exhibits a behavior with two discernible regimes: the error decreases exponentially 
for short times, then decreases polynomially for longer times. We show that the well known quadratic speedup 
over classical search is associated only with the exponential error regime. We illustrate the results through 
examples of evolution paths derived by minimization of the adiabatic error. We also discuss specific strategies 
for controlling the adiabatic error and run time. 
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I. INTRODUCTION 

In adiabatic quantum computation (AQC) OH!] quantum 
algorithms are implemented by initializing a system in an eas- 
ily prepared ground state, followed by adiabatic evolution sub- 
ject to a Hamiltonian whose final ground state represents the 
solution to a computational problem. It is known that this 
model is computationally equivalent to the standard circuit 
model of quantum computation, in the sense that each model 
can simulate the other with polynomial resource overhead [0- 
6]. While error correction methods have been proposed for 
AQC 10 Si. and arguments have been put forth that AQC is 
inherently insensitive to certain types of errors l9t- HTll . unlike 
the circuit model it is still an open question whether AQC can 
be made fault tolerant subject to realistic noise models per- 
taining to AQC in an open system setting. Indeed, even what 
constitutes a consistent picture of adiabatic evolution in open 
systems is still the subject of some debate 

01 El- Neverthe- 
less, experiments in AQC using superconducting qubits have 
made great strides recently 11411 . 

The performance of adiabatic quantum algorithms is char- 
acterized by the "adiabatic error", i.e., the fidelity loss be- 
tween the actual time-evolved state (the solution of the 
Schrodinger equation) and the instantaneous ground state, cal- 
culated at the final time T. Since the adiabatic theorem lfl5l - 
[l8ll guarantees that in the limit of arbitrarily slow evolution the 
error approaches zero, one expects that a more slowly varying 
Hamiltonian and/or a longer evolution time should result in 
higher fidelity or accuracy. This is the accuracy-time trade- 
off in quantum adiabatic algorithms. This tradeoff is often 
formalized by the adiabatic condition, which states (roughly) 
that the variation rate of the Hamiltonian should be <C the 
adiabatic error times the gap squared. A clear disadvantage of 
this condition is the inherent vagueness of "<C," which makes 
it difficult to reliably quantify the evolution time vs the de- 
sired accuracy. Moreover, violations of the traditional adia- 
batic condition have been reported OsHUlL m the sense that 
the condition is neither necessary nor sufficient for adiabatic 
evolution. While these violations have been explained as be- 
ing either due to inconsistent manipulations ll22tl or due to res- 



onant transitions ll23ll . it appears that there is no single "adi- 
abatic theorem". Rather, a number of different rigorous con- 
ditions have been derived, which apply under various math- 
ematical assumptions ll24T - t3lll . A rigorous condition which 
holds for analytic Hamiltonians and exibits the explicit scal- 
ing with system size, needed for AQC resource quantification, 
was derived in Ref . ll32tl . 

In this work we perform a comprehensive analysis of the 
adiabatic error for the case of an adiabatic quantum search al- 
gorithm, in a closed system setting. We focus on quantum 
search not only because it is an important example of a quan- 
tum speedup 13311 . but also because it is amenable to an exact 
analytical treatment. Indeed, rather than relying on a particu- 
lar form of the adiabatic theorem, our approach is based on an 
exact treatment of the underlying dynamics. We calculate the 
adiabatic error as an explicit function of the evolution time. 
We work out the formal 1/T expansion of the error — often 
used in rigorous treatments of the adiabatic theorem — through 
which we provide a large-T polynomial upper bound for the 
adiabatic error. We shall argue that a leading order truncation 
of this series expansion can result in misleading estimates for 
the scaling of the evolution time vs system size. We draw this 
conclusion on the basis of a careful study of the adiabatic er- 
ror, showing that it may exhibit a short-T exponential decay, 
which is hidden in the formal polynomial expansion. We show 
that this short-time exponential decay heralds the correct scal- 
ing for the evolution time vs system size (a quadratic speedup 
over classical search), avoiding the overestimation that results 
from the long-time polynomial decay, and which leads to a 
loss of the quantum speedup. Additionally, we propose spe- 
cific adiabatic evolution paths ("interpolations") — inspired by 
a procedure for minimization of time functionals for adiabatic 
evolution l34ll — and use these to illustrate our results. We 
also examine the applicability of the traditional measure of 
the evolution time and contrast this to what we obtain from 
our analysis. This careful analysis enables us to show explic- 
itly how one can reduce the adiabatic error as a function of the 
evolution time. 

The structure of this paper is as follows. We start with some 
general background in Sec.|II]by delineating the framework of 
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adiabatic quantum computation, the definition of the adiabatic 
error, and the adiabatic theorem. In Sec. [HI] we specialize to 
the specific problem of quantum search. There we introduce 
a general Hamiltonian interpolation for the search problem, 
and analytically solve the corresponding Schrodinger equa- 
tion. Section [TV] contains the core of our results. We first 
derive an exact expression for the adiabatic error in subsec- 
tion [IVA] Next, in subsection IIV Bl we approximate this ex- 
pression and show how the polynomial and exponential be- 
haviors emerge. Since estimation of the adiabatic error re- 
quires specific interpolation paths, we derive a general inter- 
polation for the search Hamiltonian in subsection IIV CI and 
investigate in detail three special cases. The adiabatic error 
for the general interpolation together with the special cases 
is estimated in subsection IIV D| In subsection IIV E| we dis- 
cuss a strategy for controllably reducing the adiabatic error. 
There we show how one can employ a freedom in the interpo- 
lation to manipulate the adiabatic error. This accounts for the 
performance-resource tradeoff. We conclude with a summary 
of the results and an outlook in Sec.lVl 



II. FRAMEWORK OF ADIABATIC QUANTUM 
COMPUTATION 

We begin by defining the adiabatic error as it arises in the 
context of AQC. We also provide a brief review of some perti- 
nent facts concerning the "traditional" adiabatic theorem, also 
in the context of AQC. However, we note that in the sequel we 
shall not use the adiabatic theorem; rather, we shall treat the 
dynamics directly by solving the Schrodinger equation, and 
later enforce adiabaticity by means of a 1/T expansion. 



A. Adiabatic error 

Let us treat the total evolution time T as a parameter and 
define the the scaled (dimensionless) time 

r = t/Te [o,i]. (l) 

Assume that we have an TV-dimensional quantum system that 
evolves for a total time T under the Hamiltonian Ht{t) with 
spectral decomposition 

N-l 

H t (t) = Yl ^(T)|$i(r)}(^(r)|, re [0,1] (2) 

where E.i (r) < Ej (r) for i < j, except that the ground state 
energy is separated from the rest of the spectrum by a nonva- 
nishing gap 

D(t) = Ei(t) - Eo(t) >0. (3) 

Possible level crossing among excited eigenstates 
{|$i(r))}i>o shall not concern us here because in the 
following we shall only focus on the ground state |$o(t)). 



We assume that the system is closed, i.e., the evolution is 
governed by the Schrodinger equation 

i|^r(r))=Tflr(T)|^r(r)), (4) 

where from now on dot denotes d/dr, and we set h = 1. The 
initial state is assumed to be the ground state, i.e., 

IV-t(O)) - |$ (0)). (5) 

Remark. We shall use the subscript T to indicate the para- 
metric dependence on T. When symbols already have a sub- 
script we shall avoid the additional T subscript so as not to 
clutter the notation. 

A central quantity of interest to us is the "adiabatic error" 
5^{t), which quantifies the distance between the instanta- 
neous ground state and the actual state: 

4d(r) = V1-|(1>o(t)|Vt(t))| 2 . (6) 

Note that (5 a d(r) is a distance and (in)fidelity measure in their 
rigorous sense. Indeed, the trace distance T> and fidelity T 
between two arbitrary density operators Q\ and g 2 are defined 
as Q 

T>{Qi,Qi) = ^Wei - P2II1, (7) 
F(Qi,Q2) = llv^iV^lli' W 

where ||X||i = TrVxTx is the trace norm (sum of singular 
values of X). For pure states Qi = \ipi){ipi\ it is not hard 
to show that V(g ug2 ) = y/1 - |(V>i|^ 2 >| 2 = £>(|V>i}, life)) 
and F{qi, q 2 ) = KV'ihM = -^(IV'i), l^))- Thus 

^ad(r) = V (\Mt)),\Mt))) 

= v/l-^(K'T(r)),|$o(r))) 2 . (9) 

An equivalent and useful formulation of S^t) is the fol- 
lowing. Let -P(t) = | ( l > o(' r ))( c f > o(' r )| denote the ground-state 
eigenprojection. The connection between the initial prepara- 
tion P(Q) and the time evolved state P(r) is given by 

P(r) = A(t)P(0)A1(t), (10) 

where the unitary operator A(t) — called the adiabatic inter- 
twiner i24ll — determines the evolution in the eigenspace P(t) 
and its orthogonal complement 1 — P(i~)'- 

JV-l 

A(r) = Y, |*<(r)><#<(0)|. (11) 

i=0 

One can assign a (dimensionless) "adiabatic Hamiltonian," 
Ha{t), to this evolution defined via 

iA(r) = H a (t)A(t). (12) 

On the other hand, the evolution operator generated as 

iV T {r) =TH(t)V t (t) (13) 
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dictates the actual dynamics: 

\iPt(t))(tPt{t)\ = Vt(t)P(0)V}(t). (14) 

Thus the error <5 a d ( t ) in fact measures the difference between 
Vt(t) and A(t), or equivalently how far the unitary operator 

VL t {t) = A\t)V t {t) (15) 

is from 1 (identity). Indeed, we have 

<5ad(r) = v/l-K*o(0)|fiT(r)|$ (0))| 2 , (16) 

which vanishes iff VIt{t) = 1. 

The above formulations have presumed that the degeneracy 
of the ground-state eigensubspace does not change in time. 
Nonetheless, there may be situations in which this degeneracy 
preservation assumption does not hold. In fact, as we shall 
see later, the quantum search problem we study in this paper 
falls into this category. Let us assume that the initial state is 
pure, ^^(O)) = |<£>o(0))> but the instantaneous ground-state 
eigenprojection P(r) is not necessarily rank-1, accounting for 
the possibility of degeneracy. Intuitively, if the actual state 
\iPt(t)) only has components in the support of P(r) the al- 
gorithm has achieved its goal at the instant r, whereas a less 
than full overlap denotes lack of success at this instant. The 
overlap of |^r(r)) and the support of P(r) can be quanti- 
fied, e.g., with (iPt(t)\P(t)\i/jt(t)), whence we define the 
following performance error: 



= y/1 -(Mt)\P(t)Wt(t)). 



(17) 



It is evident that < <5 ad (r) < 1, with #, d = 1 iff the ac- 
tual state has vanishing overlap with the support of the instan- 
taneous ground-state eigenprojection, while 5[ ld = iff the 
actual state resides anywhere in the support. 

We can write ($o(0)|n T (<r)|$o(0)) = Tr[P(0)n T (r)] 
and (Vt(t)|P(t)|Vt(t)) = Tt[P(t)V t {t)P(0)V^t)] = 
Tr[P(0)fi T (r)P(0)^(r)]. Using this, note the difference 
between 



and 



Mr) = V / i r PWMrf 



<^(r) = ^l-Tr[P(0)fi T (r)P(0)4(r)]. 



(18) 



(19) 



Namely, in the nondegenerate case S^t) = iff Qt(t) = 1, 
indicating that the evolution was perfectly adiabatic. In the 
degenerate case, on the other hand, algorithmic success does 
not require the evolution to be perfectly adiabatic since only 
nonvanishing overlap with the ground-state eigensubspace is 
required, i.e., 5'^ = iff V t (t)P(0)v£(t) £ supp[P(r)], 
whereas complete failure requires the dynamics to remove any 
overlap with the ground-state eigensubspace, i.e., S' id = 1 iff 

Vt(t)P(0)V}(t) £ supp[P(r)]. 

While the error as defined in Eq. ( fTTI i is not necessarily 
a distance in the strict sense, it is adequate for quantifying 
the adiabatic error. Note also that <5 ad (r) reduces to <5 a d(r) 
when the ground state is nondegenerate. From now on we 
shall be using these various expressions for the adiabatic error 
as appropriate in the rest of this paper. 



B. The adiabatic theorem 

One variant of the "traditional" adiabatic theorem lfl7ll 
states that given an < e <C 1 and a time-dependent Hamil- 
tonian H(t) with a nondegenerate ground state, the adiabatic 
error satisfies <5 ad (l) < £, provided that 



max T ||P(r)| 
min T D 2 (t) 



< eT, 



(20) 



in which j| • || is the standard operator norm, defined as the 
maximum singular value, i.e., 



|*||= sup \(v\VxVC\v)\, 

\v) ,\\v\\ —1 



(21) 



which reduces to sup^ |M|=i l( il l*l' L ')l for normal operators. 
As remarked in the Introduction, this condition is hardly quan- 
titative due to the intrinsic vagueness of "<C", and has been the 
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subject of critique (consistent with its lack of rigor) 1 19 
justifications 11221 12311 . and rigorous improvements 124143 2H 



Nevertheless, it remains a useful rule of thumb, as long as it is 
applied with appropriate care. 

An immediate implication of the adiabatic theorem is that, 
assuming it is initialized in the ground state, the system re- 
mains close to the final ground state at t = T. Thus, by 
choosing the Hamiltonian such that P(0) corresponds to a 
simple ground state |$o(0)) (simple in the sense that it is 
easily preparable), and H(l) represents a Hamiltonian whose 
ground state |$o(l)) identifies the solution to a computation- 
ally hard problem, one can devise an adiabatic version for the 
corresponding algorithmic or computational task. This is pre- 
cisely the insight that led to the advent of AQC fUd]. 

A simple "annealing schedule", or "path", between H(0) 
and H(l) is the following linear interpolation in r: 



H{t) = (1 -r)P(O) +tH(1). 



(22) 



In physical situations, however, one often realizes the dynam- 
ics by tuning some time-dependent control knobs or couplings 
x(t) = (xi(t), . . . ,xk{t)) of the Hamiltonian. This sug- 
gests that a generalization of Eq. d22l can be introduced 
by assuming access to a controllable set of non-commuting, 
linearly-independent primitive Hamiltonians {Hi}^sf N com- 
bined, e.g., as P[x(r)] = J2iLi £i (ll-^Lii, Further generaliza- 
tions can be introduced as well lul [34l l36l [3711 . However, for 
our purposes in this paper we shall consider the parametriza- 
tion HIM] 

P(x(r)) = Xi (t)H(0) + x 2 (t)H(1), (23) 
with the boundary conditions 



(*! (0)^2(0)) - (1,0), 

(x 1 {l),x 2 {l))=(0,l). 



(24) 
(25) 



Two remarks are in order regarding AQC. (i) A primary 
goal in AQC is to make <5 a d(l) decrease more rapidly for a 
given T and problem size, or alternatively, to make T smaller 
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for a given e and problem size. Often the problem size is 
given by N, the dimension of the Hilbert space. However, in 
the context of many -body quantum systems, where the Hilbert 
space is a tensor product of subsystems (e.g., qubits), log N is 
the correct measure of problem size, coinciding with system 
size, (ii) The "run time" complexity of a quantum algorithm 
should be defined as 

T im = Txmax.\\H(T)\\, (26) 

T 

not T J3]. This regularization is required because of the 
energy-time tradeoff in quantum mechanics, in the sense that 
multiplication of H{t) in Eq. d2Qb by some positive factor 
a manifests itself as dividing T by the same factor, mak- 
ing it possible to decrease T arbitrarily by choosing a suf- 
ficiently large. This tradeoff can also be understood via the 
Schrodinger equation 0, in which the final state of a system 
evolving under Hamiltonian H(t) for T is the same as that 
of a system evolving under aH(r) for T/a. This ambigu- 
ity is fixed by the definition of r lun as in Eq. d2oT ). Scaling 
of T run with system size, for a given upper bound on the error 
^ad(l) < £, determines the run time complexity of the corre- 
sponding quantum algorithm. 



where \<fr) = ^^Sq 1 \i)/y/~N is the equal superposition of 
all of the "label" states (items), P» = J2meM \ m )( m \ is 
the projection over the subspace Ai of the marked items 
= M), and J is a dimensional constant which sets the 
energy scale. In other words, the initial state \tp) is the ground 
state of the initial Hamiltonian Hi, while any state supported 
fully on A4 is a ground state of the final (oracle) Hamiltonian 
Hm ■ Note that unlike previous treatments of adiabatic quan- 
tum search J^H^H^], the Hamiltonian Hm has a degenerate 
ground eigenspace spanned by {|m)} m6 ». 

We remark that the search Hamiltonian d27l i is a member of 
the following class of projective Hamiltonians 134115(1 l5ltl : 

#(x(t))= Xl {r)Pt +x 2 {r)P( b} , (30) 

where P^ = l-|a)(a|,Pi } = l-£ b |b)(b|, with |a) and 
{|b)} fixed (normalized) vectors in the system Hilbert space, 
for which (&\P^ |a) is a given function of N. In the case of 
the search problem, we have |a) = \tp), |b) = \m), whence 
| (a|b) | = 1/ y/~N. The results of this paper can be generalized 
to other members of the class of projective Hamiltonians. 



III. QUANTUM SEARCH HAMILTONIAN 



A. Two-dimensional reduction 



Grover's quantum search algorithm 13311 performs a search 
for M "marked" items among N items of an unsorted 
database, presuming that there is an "oracle" for distin- 
guishing the marked from the unmarked items. The algo- 
rithm in its original form (M = 1) comprises the follow- 
ing steps: (i) assign orthonormal quantum states (i.e., labels) 
{|0),...,|iV — 1)} to the items, (ii) prepare the quantum 

system in the equal superposition state J^iLo 1 K)/V^> an< i 
(iii) apply the "Grover operator" — encompassing the oracle — 
repeatedly l33l l35ll . The algorithm finds a marked item after 
Trim = 0(y/ N / M) calls of the oracle — a quadratic speedup 
over the best classical algorithm — and is provably optimal for 
any N (not necessarily very large) ll38ll . Various generaliza- 
tion of the algorithm have been introduced (e.g., Refs. ll39l - 
Hl ). and it has also been implemented experimentally in a 
number of physical settings (e.g., Refs. Il45l - l49l0 . 

An adiabatic Hamiltonian version of the search algorithm 
was first introduced in Ref. 01, but failed to display the ex- 
pected quadratic speedup as it relied on the linear interpolation 
of Eq. d22l . This was fixed in Ref. ll36ll by using a non-linear 
but one-dimensional interpolation with ^(r) = 1 — x\{t) 
[recall Eq. J23H . which moves fast when away from the mini- 
mum gap, but slows down near it. As shown in Ref. (HHl, this 
result can be further improved, in the sense of a smaller adia- 
batic error, by adopting a two-dimensional interpolation as in 
Eq. (|23), 

H(t)/J = Xi (t)H x + x 2 {t)H m , (27) 

with the two projective Hamiltonians 

Hi = 1 - 10X01, (28) 
H M = 1-P M , (29) 



In the computational basis, in which \iPt(t)) 
X^q 1 V'i(' r )K)> tli e Schrodinger equation (|4]i becomes 



iipi = T 



^ m&M / j=l 



N 



(31) 

with the initial value ipi (0) = 1/ yN. It can be seen from this 
expression that the marked components all behave similarly, 
as do the unmarked components. Hence we can rewrite the 
state |-0t(t)) as 

IV't(t)) = Mr) ]T \i) + V'm(r) ^ (32) 

where the subscripts "u" and "m" denote "unmarked" and 
"marked," respectively. The normalization condition now 
reads 

(N - M)|^ u (r)| 2 + M\M{t)\ 2 = 1. (33) 
By defining the (unnormalized) two-dimensional vector 

= (^u>m) T , (34) 

and the (non-Hermitian) reduced Hamiltonian matrix 



Hj J 



in which 



rx\ + x 2 —rx\ 
(r - 1):ei (1 - r)xi 



r = M/N 



(35) 



(36) 



5 



is the fraction of the marked items, the Schrodinger equation 
reduces to 



i\Mr))=TH{T)\Mr)), 



(37) 



with the initial condition \tp(0)} = (1/VN, l/VN). There- 
fore, not only is the parameter space of the problem two- 
dimensional, it is described by an effectively two-dimensional 
Hamiltonian (in the m-u representation). This reduction from 
the real Hamiltonian H [Eq. (I27H to the effective Hamiltonian 
H [Eq. (l35l ll will prove useful in our analysis below. 

Later in the paper we shall need the norm of the Hamilto- 
nian as well. This can be calculated easily from Eqs. (|2TT > and 
}, and yields: 



1^)11 = 



and similarly, 



\H{r) 



J\ xi{t) + x 2 (t)\ : 
JsJT^ T\x x {t)\ : 



J\±i(t) +x 2 (t)\ : 
Jx/T - 7|ii(r)| : 



Xi(t) + x 2 (t) t^O, 
otherwise, 



(38) 



±i(t)+x 2 (t) ^ 0, 
otherwise. 



(39) 



B. Diagonalization and unitary interpolation 



The Hamiltonian H/ J d27l i has three distinct dimensionless 
eigenvalues E- < E + < E> , where 



where 



E T = (xi + x 2 T A)/2, 
E> = xi + x 2 , 



A = sj (xi — x 2 ) 2 + Arx\X2, 



(40) 
(41) 

(42) 



is the dimensionless gap (hence D = J A) and E > is (N — 2)- 
fold degenerate. 

Let (T z = diag(l, —1) and a y = (° denote the Pauli 
matrices. Let us define the similarity matrix 

S 

and the unitary 



_ ( r/y /l-r ^/r 

-VI —r y/f J ' 



(43) 



— p—ivy arccos[(xi — (1 — 2r)x2)/A]/2 



(44) 



Then the effective Hamiltonian H(t) [Eq. d35l ll satisfies 

H/J = SH(t)S'\ (45) 

where 

H(t) = ldiag(-B+, E^)A^ (46) 
is the Hermitian core of H / J, and we easily find that 

H(t) = i {[xx(r) +x 2 (j)]l + A(r)A(r)a 2 lt( T )| . 

(47) 



This last result is remarkable: it states that, up to an overall 
(time-dependent) shift [x\ (t) + x 2 (t)] 1 and a conformal fac- 
tor A(r), the reduced Hamiltonian H{t) is a unitary inter- 
polation A(r)a z A^ (t) ll52tl . We will exploit this observation 
below. 

The non-Hermitian reduced Hamiltonian H(t) and its Her- 
mitian core H{t) have the same set of eigenvalues E T = E T , 
and we have the spectral resolution 



where 



H/J = E- !*_)(*_ I + (* H 



\<$> t {t))=A{t)\z,t), 



(48) 



(49) 



and |$ip(0)) = \z, =f) are the eigenvectors of cr z , correspond- 
ing to the eigenvalues =pl. The unitary operator A{t) acts as 
a reduced adiabatic intertwiner [Eq. (TTOlil. in the sense that for 
the reduced projection P T (r) = |$ =f (t))($ =f (t)| we have 



P T (r) = A(r)P T (0)^(r) 



(50) 



Remark. We emphasize that throughout the paper hat and 
tilde denote states or operators in the reduced representation; 
the only exception is A. 



C. Solving the Schrodinger equation 

In solving the Schrodinger equation and calculating the adi- 
abatic error 5 a d(l), it is more convenient to work with the nor- 
malized state 

I^^JeVmS-^tM). (51) 
Equation (|37T i now becomes 

i$ T (T))=TH( T )\MT)), (52) 
with \i/jt(0)) = \z, — ). Solving this equation results in 

\Mr)) =V t (t)\M0)), (53) 

in which 



V t (t) = Texp -iT I H(t') dr' 



(54) 



is the time-ordered reduced evolution operator. 



1. General setup: Adiabatic interaction picture and Dyson series 

Having observed that the Grover search problem can be cast 
as a conformal unitary interpolation, we outline a general, 
systematic approach for solving the corresponding class of 
Schrodinger equations using the Dyson series expansion ll53ll . 
Consider as a specialization of the general time-dependent 
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Hamiltonian of Eq. (O the "conformal unitary interpolation 
Hamiltonian" 



H t (t) = A(t)A(t)H t (0)A^i 



in which 



N-l 



(55) 



(56) 



i=0 



is the spectral decomposition of the (traceless) initial Hamil- 
tonian H(0), the unitary operator A(t) satisfies the adiabatic 
Schrodinger equation (fT2] > generated by the adiabatic Hamil- 
tonian H a (t) Il24[l50ll . and (the dimensionless gap) A(r) > 
is a smooth function with the initial value A(0) = 1. It is 
also useful to think of the time-dependent Hamiltonian Ht{t) 
[Eq. ( TSBI ll as the "adiabatic interaction picture Hamiltonian", 
though normally an interaction picture transformation does 
not involve a time-dependent prefactor such as A(r). It is 
evident that the eigenvalues and eigenprojections of Ht{t) 
satisfy 



Ek(r) = A(r)£ 4 (0), 
P,(r) = A(T)Pi(0)A\T). 



(57) 
(58) 



Equation d57l i implies that the spectrum of Ht{t) evolves 
conformally with A(r), i.e., all the eigenvalues are multiplied 
by the same factor, while Eq. ( l58l implies that the eigenpro- 
jections are unitarily connected, and the degeneracy is con- 
stant in time. Comparison of Eqs. ( TTOb and d58l reveals that 
P(t) = Pq(t) and here too, A(t) plays the role of the adia- 
batic intertwiner. 

Rather than solving the Schrodinger equation 



Vt(t)=TH t (t)V t (t), 



(59) 



we solve the equation of motion for the "adiabatic interaction 
picture propagator" [cf. Eq. ( fTBI ll 



Qt(t) = A\t)V t {t). 



(60) 



This provides a more direct tool for the calculation of the adi- 
abatic error [Eq. ((6)]. Let us define 

H (t) = A\t)H t {t)A(t) = A(r)ff T (0), (61) 
Hi(t) = A^(t)H a (t)A(t). (62) 

Note that Ho (r) has dimensions of energy while Hi (r) is di- 
mensionless. In the adiabatic interaction picture the (dimen- 
sionless) "perturbation" is THq(t) — H\(t), i.e., it follows 
from Eqs. (O and (|59]l that 



itl T ( T ) = [TH (t) - Hi(t)]SIt(t) 



(63) 



We also define the two unitaries Vq (t) and V\ (r) through the 
following equations: 



iV = TH V , 



(64) 
(65) 



It is easily seen that VqVi also satisfies Eq. (l63l . so that 

Qt(t) = V (t)V 1 (t). (66) 

To simplify the analysis, we only consider Hamiltonians for 
which 



H a (t) = h A (r) S, 



(67) 



in which h A is an integrable function and H is a constant (r- 
independent) and traceless operator belonging to the space of 
linear operators acting on the system Hilbert space. Thus from 
Eq. (l62l l we obtain 

(68) 



Hi = H A . 
Note that, from Eqs. d6Tb and ( |64l i, 

V (r)=e- lTH ^r A , 



(69) 



where J T A is shorthand for J Q r A(t') dr' — we shall use the 
similar shorthand 

J g = J g(T f )dT f (70) 

wherever convenient. Inserting Vq into Eq. ( |65] l yields 

iV^/CrVi, (71) 
in which the kernel JCt is defined as 



K, t {t) = ih A (r) e 



\TH T (0) J T A n -iTH T (0) J r A 



Equation (fTTT i. or equivalently the Volterra equation 



(72) 



(73) 



Vi(t) = 1+ / K, t {t')Vi{t') dr' , 
Jo 

can be solved iteratively, yielding the Dyson series 

OO /»T /'Ti-l 

V^(r) = 1+5^/ ACr(n)dn... / ACr(77)d77. 
i=1 Jo Jo 



(74) 



2. Quantum search 



Now we apply the method described above to the adiabatic 
quantum search problem. Recall that we are working with the 
reduced states and operators (hence the hat over all reduced 
quantities, except the gap A [Eq. (1421 1). Comparing Eq. 
with Eq. (EBT l implies that 



H(0)/J= 2 CTz 



(75) 



after excluding the trivial term oc 1 from H. Similarly, 
Eqs. gU and (|67]) yield 



el 

d7 



■ arccos[(a;i — (1 — 2r)x2)/A] 



03 



XiX 2 - XiXi 

A 5 ■ 



(76) 
(77) 
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As a result, from Eqs. d69l and (l72l we obtain 



% = e-«diag(e-i- /T J' rA > e* jr J' TA ) > (78) 

(79) 



where the phase factor 



X(r) = \ JT l [*i{t')+x 2 {t')] At' (80) 

compensates for the removal of the trivial term from Ht{0). 
Some simple algebra then yields 

where, for £ > 1, the Dyson series terms are 

2|(r)= r^Cr'^-iCrOe^-^'^^^dr', (82) 
Jo 

and Io( T ) = 1- This completes the derivation of 

n T =A^V = %Vi. (83) 

IV. ADIABATIC ERROR IN THE SEARCH ALGORITHM 

In the previous section, we worked out the solution to the 
Schrodinger equation in the quantum search problem. Having 
collected the pertinent ingredients, we now return to calculat- 
ing our main object of interest, <5 a d(V). 

After deriving an exact formula for the error, we proceed 
with approximating it in the large system-size limit, identified 
with r <C 1. We start with the well-known polynomial ex- 
pansion of (5 a d(l) in terms of T, which works well for large 
times. In refining this result, we show that in fact two regimes 
are discernible in the behavior of (5 a d(l) vs T: (i) the onset of 
exponential decrease, followed by (ii) a polynomial tail. This 
dichotomy will appear to be crucial in a correct characteriza- 
tion of the scaling of the run time of the algorithm with system 
size. 

The exact behavior of the error in the algorithm depends 
strongly on the form of the interpolation one chooses for the 
Hamiltonian. Inspired by our earlier study aiming at minimiz- 
ing the adiabatic error in quantum algorithms IdOL we shall 
suggest a general class of interpolations, which includes three 
specific cases already studied in the literature. Next, we in- 
vestigate the specific behavior of the adiabatic error for each 
interpolation, separately. Finally, we shall suggest methods 
for suppressing the adiabatic error even further. 

A. Exact relation 

Recall that the adiabatic quantum search Hamiltonian (|27T i 
has a nondegenerate ground state \(f>) at the initial time r = 0, 



whereas the ground-state eigenprojection at the final time r = 
1 is Pm-i which is 7A/-fold degenerate. Here, any full superpo- 
sition of the form ^^(l)) = J^ieM V'iK) will work equally 
well, whereas if does not have complete support over 

Pm then this indicates that the algorithm has partially failed. 
Hence, following the discussion in Sec.[II][Eq. (fTTIil, the adi- 
abatic error at the final time is determined by 



Mi) = Vi-(^)\Pm \Mi)} 
Vi- W m (i)| 2 

ED 



aVT^#u(i)I- 



(84) 



An equivalent formulation can be obtained for the two- 
dimensional reduction we discussed in Sec. [HI] In this rep- 
resentation |V>r(l)) is replaced by |V>t(1)) = V^(1)|$_(0)) 
[Eq. (|53ll; similarly, the instantaneous ground state is repre- 
sented by the nondegenerate state |$_(1) = A(1)|<1>_ (0)) 
[Eq. (|49ll. Thus we can employ the error formula appropriate 
for the non-degenerate case [Eq. (fISVl. whereby 



Ml) 



i-|<<Mi)IV>t(i)>P 



i-\( z ,-\n T (i)\ z ,-)\ 2 
K«,+|nT(i)K-)| = K*+(i)|$r(i)>|,(85) 



in which S1t(1) is given by Eq. (1831 ), and in the last line we 
used the unitarity of Q,t(1). The equality of Eqs. d84l l and 
is immediately seen from 



K#+(l)$r(l)>| 1> 

{43},g4} 



M 



(io)lt (1 ) 5 ^(j : «) 



JVvT^7|^u(i)|. 



(86) 



Inserting S1t(1) — noting Eqs. d78l ) and d8lT l — into Eq. ( T85l > 
yields the following exact expression for the adiabatic error at 
the final time: 



mi) = 

which is upperbounded by 

oo 

Mi) <E l^+i (i)|- 



(87) 



(88) 



From the above equations we can in principle calculate the 
adiabatic error or its bound given that we know all I Q dd(l)'s. 
Remark. The value of (5 a d(T = 0) will be important later. 

In this case, we have |V>t=o(1)) = thence Eq. 

yields 



MO) = V^r, 



(89) 



in which we used the boundary conditions ( fZ5b in A(l). This 
relation is valid for any interpolation that satisfies the bound- 
ary conditions. 
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B. Approximation of the adiabatic error 



obtain 



An exact calculation of the adiabatic error from Eq. (l87T i 
can be challenging because of the infinite number of terms in 
the summand and the fact that each term contains a multiple 
integral. To alleviate this difficulty, in this subsection we ap- 
proximate the upper bound on <5 a d(l) from Eq. (|88l from the 
first few X;'s, and argue that this suffices for most algorith- 
mic purposes. We start from an expansion in powers of 1/T, 
based on integration by parts, and explain its limitations. We 
then provide more careful analyses, based on the residue the- 
orem and on the stationary phase method, both of which lead 
to an exponential error estimate. 



1. Polynomial expansion 

The most common rigorous adiabatic approximation em- 
ploys an ex pan sion in powers of 1/T, presuming that T is 
"large" 1I24I l29l l32tl . Let us now show how one can system- 
atically expand the adiabatic error as a polynomial inl/T by 
extracting powers of 1/T from each term Ii (t) through inte- 
gration by parts. 

From the identity e~ lTY ^ = i/ (TY(t)) £e~ lTY ^ 
[valid for any differentiable function Y(t)], we obtain the fol- 
lowing relation by integration by parts: 



o 

i 

TV 



G (r') e - iTY (r') dT ' 



(90) 



Y(r'} 



dr' \Y(t') 



Notice how this extracted a 1/T in front of the first term. 
In the second integral on the right hand side we can iterate 
the same trick of replacing the exponential with its derivative; 
which generates 1/T 2 and higher order terms. This provides 
a systematic way for generating poly(l/T) expansions of ex- 
ponential integrals, as we shall see more specifically below for 
the Z/(t)'s. 

Using Eq. d90l . we obtain 



Ii(r) = — 
; JT 



-iJTj T ' a M t ') 

A(r') 

d 



-iJT j T A 



(91) 



Applying once more the exponential identity < f90l > for the sec- 
ond integral above gives rise to 0(1/T 2 ) terms; whence, 



Mr) = — 



h A (r) 



iJT } T A 



A(r) 



MO) 



(92) 



i.e., |Ti(t)| = 0(1/T). In analogous fashion, for I2 (r) we 



Mr) = ± 



T h\(r') A _, 



A pfdr'-M0)2?(r) 

T dT'MrV JTr ' A 
T ' dr „ d ( h A {r") \ r ^., T r' 



dr" V A(r") 



(93) 



from which 



i 




h A (0) r 


If 


-Jo A(r') . 


(JT) 2 l 



A(r) 



-h A (0) 



O 



(t 3 )' 



(94) 



and |T2(t)| = 0(1/T). By induction, one can conclude from 
Eq. <H3) that Sill 

1 \ 



\X2i-i(r)\ = \M(r)\=0^ ¥ ) 



for / e N. Thus, from Eq. 
becomes 



(95) 

the adiabatic error bound 



Ml) < |2Ti(l 

1 r 

JT 



1 



0( T^) 
|^(0)| + |^(1)| 



(96) 



This relation can be simplified further. From the bound- 
ary conditions (|24]i and (f25]l, we obtain A(0) = A(l) 
i [ Eq. 63] , h A (0) = y/r(l-r)x 2 (0) and h A (l) = 
-^r{l - r)ii(l) [Eq. d76tl. Thus Eq. d%]l reduces to 



<5ad(l) < 



JT 



O 



1 

T 2 



(97) 



Proceeding in a similar manner, one can in principle obtain the 
exact form of the coefficient of each 1/T' term, for arbitrary 
I G N. 

Remarks. Let us make some remarks regarding the polyno- 
mial expansion, and in particular Eq. d97l ). 

(i) Notice that often the poly (1/T) series is truncated af- 
ter the first or at most the second order term, on the ba- 
sis of the assumption that for sufficiently large T the first 
couple of terms should give a reliable and accurate upper 
bound. However, without correctly defining what "large" T 
means, a truncation after the first few terms might be unjus- 
tifiable. In fact, in addition to T, the system size log N (in- 
troduced here through r) and the gap A(r) are also key play- 
ers in the estimation of S^t). The minimum (system-size 
dependent) gap A m j n = min T A(r) works in general as a 
bottleneck for t he p erformance of quantum algorithms (e.g., 
Refs. |H HHl HI IH). Specifically, where the gap closes 
or becomes small, the adiabatic approximation may not hold, 
indicative of a "q uantum phase transition" (in the thermody- 
namic limit) l56U59ll . This implies that the coefficients 
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of some high order 1/T 1 terms might have a stronger gap de- 
pendence than those of lower order terms. With this caveat, 
neglecting those higher order terms is not always possible. In 
fact, it is not difficult to see that the coefficient of the 1/T 2 
has a A~ 6 dependence (see also Ref. I29I0 . stronger than the 
A~ 2 dependence of the coefficient of the 1/T term in S^t) 
[Eq. 476)]. 

(ii) For similar reasons, an estimate of T arising from 
<5ad(l) < £ (for a given e) along with a truncated poly(l/T) 
expansion, is not always reliable. We shall see this explicitly 
later in this section. 

(iii) One might argue that <5 a d(l) = 0(1/T) results from 
an energy-time uncertainty relation such as <5 a d(l) x T ps 1 
[presuming <5 a d ( 1 ) is directly related to the uncertainty in mea- 
surement of energy]. However, this argument is not rigorous 
and should not be considered as a replacement for the analy- 
sis leading to the poly(l/T) expansion (unless justified care- 
fully). A rigorous energy-time uncertainty relation is given, 
for example, by the Mandelstam-Tamm inequality 



A^[H] x T^K] > 1/2, 



(98) 



in which T^K] = A 4 ,{K]/\d{ip\K\^) /dt\, K is any observ- 



able, and A^\X] = ^J{ip\X 2 \il)) - {^\X\tP) 2 |6fl]. Hence, 
as is well known the naive relation A.^ [H] xTkI should be 
treated with care. Strictly, a relation between /jj.ff(l)] 

and <5 a d(l) can be constructed as the following. Note that we 
can write 



\Mr)) = y/1 ~ Mt)|S_(t)) + WI^W), (99) 

where |$^(r)) is orthogonal to |$_(t)) [Eq. |85]. Hence, 
after some algebra we obtain 

^(^(1) « <5ad(l) 

2E-[E-($±\H*\Q±) @±\H\$±)]\ r=1 + 0(51). 

(100) 

Despite this relation, connecting T and T^ t ^[K] is not 
straightforward. Although usin g di fferent versions of the 
energy-time uncertainty relation l!6lT - l63ll may provide addi- 
tional insights, we shall not further pursue this here. 



2. Exponential estimate 

a. Residue theorem analysis In the previous subsection 
we used integration by parts to arrive at a polynomial ex- 
pansion. Let us now show that an alternative, more care- 
ful analysis based on the residue theorem of complex anal- 
ysis, reveals that the adiabatic error decays exponentially for 
sufficiently short times. In some sense, this exponential be- 
havior is reminiscent of the well-known Landau-Zener for- 
mula for two-state quantum systems l64ll . which — in its sim- 
plest form — states that the tunneling probability Pt(1) = 
|^(0)|l^(l)|*i(l))| a from |^(0)) =|*o(0)} to ^(l)), for 
the Hamiltonian H(t)J J — (rer z + Da x )/2, is 



Pt(1) 



-irJTD 2 /2 



We notice that pr(X) is in fact intimately related to 5 a d(l) in 
this simple two-state case, as 



PT (l) = l-5 2 d (l). 



(102) 



The exponentiality of the adiabatic error vs time and (some 
power of) the gap has been previously established ingeneral- 
ity in rigorous treatments of the adiabatic theorem ll26l - l28l[32ll 
(see also Ref. lf55ll in a more restricted setting). 

In what follows we focus on T\ (t); higher order terms can 
be treated similarly. We start from Eq. d88l , whence 



Mi) < 



h A (r)e UT ^ A dr 



(103) 



Here ". . . " denotes the higher order terms |22i+i(l)| (I > 1), 
whose neglect we justify below in the specific examples we 
discuss. An exponential error term can be obtained, for ex- 
ample, by extending the integral to the complex time plane 
and using an appropriate closed contour for the application of 
the residue theorem l65ll . A precursor to this treatment of the 
adiabatic error can be found, e.g., in Ref. lf55ll . 

If h A (z) (z € C) is not constant it has poles at points zq 
where the gap vanishes: A(zo) = [Eq. (|76ll1. From Eq. d42l . 
we obtain xx(zo) = £2(^0) = 0, or (assuming x%(zo) 7^ 0) 



X21 



Oo±) = l-2r±2iy/r{l-r), 



(104) 



where X21 = X2/X1, The poles zq can in principle be obtained 
by inverting this relation for a given interpolation £21 (t). 
Note, however, that there might exist other singularities (typ- 
ically at infinity) arising from the exponential e~ lJT 'S A in 
Eq. ( 1103b . Therefore, estimating the integral ( 1103b requires 
finding all contributing singularities in a suitably chosen con- 
tour C in the complex r-plane. 

The value of the integral | J Q h A (T)e %JT $ A dr| can now 
be obtained by calculating the residues of the integrand at the 
poles enclosed inside the contour, 



2iti Res 

2o£inside(C) 



h A U 



JJTf'A 



C-[0,1] 



h A {z)e iJT A dz 



(105) 



For example, if the poles z$ are simple and arise from 
A(zq) = 0, calculating the integral becomes straightforward. 
Recall that h A oc 1/A 2 [Eq. ( |76T >1 and the residue of a ra- 
tional function P(z)/Q{z) at a simple pole z* is given by 
P(z*)/d z Q(z*) [65]. Thus, in this case we obtain 



Res 



h A (z)e 



UTJ'A 



•-'0 



(106) 



VKr^)x 21 (zo)e' JT T° A 
±[{l-x 21 (z)) 2 +irx 21 (z)}\ 



]_JJTf*° A 
4/ 



(101) 



where in the last line we assumed ±21(^0) 7^ 0. Note that this 
computation of the residues does not necessarily hold when 
h A = const., or when zq is at infinity. 
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This yields that an exponentially decreasing contribu- 
tion emerges from Im[J 2 ° A] of the exponential within 
the integrand — perhaps in addition to a generically non- 
exponentially-decreasing term. Thus, from Eq. (1106b we find 



Mi) 



< 



2o£inside(C) 



Rt(zo) 



■JT lm[j"° A] 



#t(1), (107) 



where Rt(zo) G K + results from the non-exponential contri- 
bution of the residue at zq and R' T (1) G K + encapsulates the 
collective nonvanishing contribution of other segments of the 
contour as well as other non-exponentially-decreasing contri- 
butions emerging from the ". . . " terms in Eq. ( 11031 ). 1 

b. Stationary phase analysis An alternative and comple- 
mentary approach for obtaining the exponential contribution 
to the adiabatic error is to use the stationary phase method. 
This method is appropriate for obtaining asymptotic value 
of complex integrals of the form J F(z)e lTG ( z >dz, where 
T > is a large number and 7 is a path for the integration. 
Assuming F(z) is a slowly-varying function over 7 and G(z) 
is an analytic function, one can argue that the major contribu- 
tion to the integral comes from the point(s) zq at which G(z) 
has a minimum, whence |6q] 



This is a less stringent condition than the standard condi- 
tion < f2Qb or its more rigorous counterparts |24M32ll . involving 
higher powers of the gap. The crossover point T* between the 
exponential and the polynomial regimes can be estimated by 
solving 

JMr) e-^M/- A] „ ^(^(1)1 + ^(o)!), (HI) 

in which Rr(r) is a non-exponential prefactor given by 
Eq. (fT07T> or dI09l . If T < T* (T > T*) the exponential 
(polynomial) behavior prevails. 

Later in this section, we show explicitly that the expected 
run time for the quantum search algorithm is often given 
by the exponential contribution; the polynomial regime may 
overestimate the minimum run time required for reasonable 
accuracy. 

Since the scaling of the run time depends on which in- 
terpolation we choose for the Hamiltonian, in the following 
we shall obtain specific interpolations by employing some re- 
cently developed results for (partial) minimization of the adi- 
abatic error IMS. 



F( z y TG ^dz 



2 71 



TG(z ) 



F(z )e 



in/4 e iTG(z ) _ 



(108) 

To apply this result to Eq. (11031 ), we replace G(z) and 
F(z) with J J z A and 1ia(z), respectively. Accordingly, zo 
is where G{z<yj = A(zo) — 0. The minimum point zq is often 
a complex number. In fact, in quantum many-body situations, 
the gap A is a nonnegative function, often with a nonvanishing 
minimum, which becomes zero in the thermodynamic limit 
J56ll . In such cases, G(z) does indeed have a minimum. It is 
also required that h a{zo) be finite [this of course is not sat- 
isfied when Ha has a pole or singularity at Zoi- If all these 
conditions are satisfied, Eqs. ( 1 1 031 > and (IIO8I 1 yield 



Mi) £ 



ftytOo)< 



2- 



JTA(zq) 



,-JTIm[/*o A] 



(109) 



as we wished. This relation complements Eq. d 1 06b in that it 
may be applicable when Eq. ( 1106b is not. 

c. Discussion The emergence of the exponential depen- 
dence of <S a d on T [Eqs. ( 1 1071 ) and (11091 )1 in contrast to the 
polynomial dependence [Eq. d97| )l is remarkable, as it indi- 
cates a much faster decay of the adiabatic error than what is 
suggested by the standard 1 /T expansion. In the exponential 
regime it suffices that T be large compared to 1/Im[/ Z ° A], 
or roughly i55ll : 



1 

JA' 



1 In the case of Eq. )106t . and from Eq. fll05t . we obtain: Rt 
is in good agreement with Eq. (89) when T = and r«l. 



(HO) 



tt/2. This 



C. Hamiltonian interpolation 

The set of available control knobs (x), as well as the way 
one varies them, determine the specific Hamiltonian inter- 
polation implemented in a laboratory setting. Theoretically, 
though, there are various ways by which one can obtain fami- 
lies of Hamiltonians for an adiabatic quantum algorithm. One 
natural choice is interpolations which minimize "physical" 
cost. For example, in the setting of Refs. ll34l[36ll . time func- 
tional were constructed from a local version of the adiabatic 
condition (1201) , whose minimization resulted in a set of Euler- 
Lagrange equations for the underlying interpolations. A dif- 
ferent method was suggested in Ref. [503, where it was shown 
that in the standard poly(l/T) expansion of the adiabatic er- 
ror 6 3 a(t) [general counterpart of Eq. d9"6l > or d9~71)l, the coef- 
ficient of the 1/T term has a geometric part, in a differential 
geometric sense. Specifically, this geometric coefficient, in 
terms of the ground-state eigenprojection, P(r) is 



[P(x(r')),P(x(r'))]||dr'. 



(112) 



Minimizing this coefficient yields adiabatic "geodesic" inter- 
polations, which partially decrease <5 a d(i~) for a given T. It 
follows from standard variational calculus l65l l67ll that the 
geodesic interpolations satisfy the following equation 



|[P(x(r)),P(x(r))]| 



o(r) 



= const. 



(113) 



where the constant is chosen to satisfy boundary conditions. 
We adopt this method in the following and derive geodesic 
interpolations for the adiabatic quantum search. 
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1. General case 

In the reduced two-dimensional representation, P is re- 
placed with P_ = |$_)($_|. Thus we have 

[P-,P-\ = + 

x(($_|$_) - ($_|$_)). (114) 
Noting with the help of Eqs. d44l . (|49T >, and (|76]l that 

(^.\$.) =ih A (z,-\A^a y A\z,-) = 0, (115) 

we find 



[P-,P-}\\ = V2<*-|*-> 



(116) 



Hence, according to Eq. dl 13b . a class of adiabatic geodesies 
can be obtained from 



/i^(x(r)j = const. = ip, 



(117) 



which implies that the adiabatic Hamiltonian H a (t) is in fact 
constant. In other words, from A(t) = e~ l f Ha [Eq. (l44l l 
we can see that 



A(r) 



-IXfTOy 



(118) 



This equation suggests that a more general choice for the in- 
tertwiner A can be obtained by (pr — > (p9(r), for some ar- 
bitrary differentiable 9; i.e., A(r) = e ~ l ^ e ( T )' T ». In terms 
of h A , this translates into choosing h A (x(r)) = ip9(r). In 
terms of the geometric factor ( II 121 ). this simply means that on 
the geodesic, the time is determined by 9(t) rather than r. To 
see this, note that 



\[P{t%P{t')]\\&t' = 



0(r) 



\[d e P(9),P(9)]\\ d6, 
(119) 



which in turn, from variational calculus, gives the following 
equation for the geodesies: 



\\[d e P(x(9(T))),P(x(9(T)))} 
Hence, 

\\[d e P-,P-]\\ = 



x gcu (e(r)) 



const. (120) 



j||^->-P-]| 



e 

const., 



(121) 



where we assumed 9 (t) > 0, 9(0) = 0, and 9(1) = 1. As a 
result, we choose the adiabatic geodesic satisfying 



(122) 



h A x(0(t)) = <p6(T) 



We remark that introducing an arbitrary nondecreasing 
function 9(t) will serve as more than a generalization. In fact, 
we show below that 9 allows us to substantially enhance the 
suppression of the adiabatic error. For example, choosing a 9 
such that it is a differentiable function of r (up to some con- 
trollable order, say, k), with vanishing derivatives (up to the 
same order k) at the initial and final times, can suppress the 
adiabatic error [up to poly(l/T fe+1 )]. This property — which 
comes at the relatively small price of sufficient control at the 
beginning and end of the dynamics — may have immediate ap- 
plications in experimental realizations of quantum annealing 
and AQC. 

The above geodesic equation can be solved analytically. 
Recall that X21 = X2/X1, with the boundary values X2i(0) = 
and £21 (1) = 00 [Eqs. (GUi and dBVl. Then from Eq. dT6> 
we find 



y/r(l - r)±2l(T) 



[1 -x 2 i(t)} 2 +4ra; 2 i(r) 
The solution to this equation can be written as follows: 

sin[29(T)(p} 



X2l(9(T)) 



sin[2(l - 9(r))(p] ' 



where we must choose 



ip = arctan[i/ (1 — r)/r\ . 
Observe that X21 (9(t)) has the following symmetry: 

o; 2 i(1-0(t)) - 1/xai (0(r)), 
which is satisfied, for example, by requiring 
x 2 (6(t))=x 1 (1-6(t)). 



(123) 



(124) 



(125) 



(126) 



(127) 



Equations ( 1124b and j 127b identify a two-dimensional inter- 
polation for the quantum search Hamiltonian. 

Notice that, given Eq. ( 11241 i. we can also add another re- 
lation (satisfying the boundary conditions) so as to find other 
classes of interpolation. For example, we can choose 



Xi (9(t)) + x 2 (9(t)) =o(t), 



(128) 



in which a(r) can be a smooth function with the boundary 
values a(0) = a(l) = 1; e.g., a(r) = 1 + r(l — r) or 
a(r) = 1 + sm(7TT). Choosing a form for a(r) corresponds 
to assuming a given ||i/(r)|| [Eq. (l38ll. which implies a given 
control over the maximum amount of the available energy in 
the system. Thus, Eqs. (1124b and ( 1128b yield 



Xl (9(r)) 



x 2 (8(t)) 



a(r) sin[2(l - 0(r))v?] 
2sin(» cos[(l - 29(r))ip] 
a(r) 



1 + cot(yj) tan[(l - 29(r))tp) 
a(r) sin[20(r)<p] 



2sin(<^) cos[(l - 29(r))<p] 
a(r) 



1 - cot(v3)tan[(l - 29(r))ip] 



, (129) 



. (130) 
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Note that this interpolation also satisfies the symmetry (1127b . 

It is clear that one can consider other auxiliary or con- 
trol conditions over the Hamiltonian different from Eq. (1128b . 
In the following, we address three special cases: (i) The 
Hamiltonian interpolation is linear in time, namely, x\{t) = 
1 — x 2 (t) = 1 — r; (ii) Hamiltonians with constant norm; 
specifically, \\H(t)\\/J = 1; and, (iii) Hamiltonians with con- 
stant gap; specifically, A(t) = 1. 



2. Linear interpolation 



If we choose 



(r) = \ - — arctan[(l - 2r)tan^l, (131) 
2 zip 



and assume a(r) = 1, from Eqs. ( 1129b and J 1 30b we obtain 
the simple linear interpolation 

Xi {t) = 1-t, (132) 
x 2 {t)=t. (133) 

3. Constant-norm interpolation 

Let us assume < Xi,x 2 < 1. The constraint 
\\ H ( T )\\/J = 1 implies that x\ + x 2 = 1 [Eq. d38ll. or 
equivalently a(r) = 1 [Eq. ( 11281 )1. Hence, in this case from 
Eqs. ( 11291 ) and ( 1130b we obtain 



x 1 (9(r)) = 
*a(0(r)) = 



1 



2 2VT^ 
1 Jr 



2 2y/T^r 



tan[(l-20(r)) V ], (134) 
tan[(l - 20(t))lp\. (135) 



This interpolation is a gene ralization of the interpolation ob- 
ill by using a local adiabatic condi- 



tained in Refs 
tion. 



4. Constant-gap interpolation 



Rather than assuming condition ( 1128b . here we consider the 
case in which the gap is constant, e.g., A(r) = D(t)/J = 1. 
Hence combining the following relation [Eq. (I42b l: 

Xl {e{T)) -x 2 {6{t))^ + Arx 1 {B{T))x 2 (0{T)) =1, (136) 
and Eq. ( 1124b yields 

Xi (6(t)) = J. sin[2(l-fl(r))y], (137) 

2y/r{l - r) 

x 2 (8(t)) = J Bin[2fl(T)y]. (138) 

2 v /r(l - r) 



D. Error estimation for different interpolations 

Having given a general recipe for adiabatic interpolations 
and having obtained three particular interpolations, we pro- 
ceed to compute the adiabatic error for each of these interpo- 
lations. Our main interest here is to analyze how the run time 
scales with system size for each of these three interpolations. 
We shall also contrast the result for each case with the estimate 
obtained from the traditional adiabatic theorem. As a result, 
we will demonstrate that the traditional adiabatic condition is 
not always reliable for estimation of the minimum run time 
(given an error threshold), whereas the time we obtain from 
the exponential regime of the adiabatic error is in fact accu- 
rate. A remarkable feature of this result is that the estimated 
times (in an appropriate sense) need not be very large. 



1. Linear interpolation 



In this case, (xi(r), x 2 (r)) = (1 — r, t), so that we have 



\\H(r)\\/J 1 1, 
\\H(r)\\/J 1 1, 



A, 



G3 



(139) 
(140) 
(141) 



Hence the traditional adiabatic condition d20b implies that, for 
i5ad(l) < s, we should have 



Jer ' 



or equivalently from Eq. ((26]), 



= - 



(142) 



(143) 



That is, that the adiabatic quantum search with a linear in- 
terpolation Hamiltonian requires a run time O(N) [recall 




1000 



FIG. 1: (Color online) 5^(1) for the linear interpolation x(r) = 
(1 — t, r), obtained by numerically solving the corresponding 
Schrodinger equation, for r = 0.05. The dashed lines represent 
the exponential fits e~ 7TjTr ^ 8 [Eq. J147H and the polynomial fit 
2y/r/{JT) [Eq. J144H . from left to right, respectively. 
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r = M /ATI, hence performing no better than a classical search 
algorithm 1 3^] . 

Note, however, that if we truncate the poly(l/T) expansion 
( |97l i after the first term, in the r«l limit we obtain 



Je ' 



(144) 



i.e., T run = 0(yN), which is not the right estimate. This 
illustrates the caveat we discussed in Sec. lIVI 

Now we employ the results we developed in Sec. [TV] for 
estimating the exponential regime of <5 a d ( 1 ) [Eq. (I107H . It is 
obvious that for this case there exist no point at which both 
Xx(t) = 1 — r and — r vanish. Thus, the pole zq (where 
the gap vanishes) is obtained simply by inverting Eq. (11041 ). 
i.e., 



1 

= - ± 
2 



2y/T^T' 

The integral over the gap can be evaluated explicitly as 



Im 



A(z) dz 



(145) 



(146) 



For the integral ( 11061 ). we choose the contour C to be a rect- 
angle composed of: (i) C\ the real line [0, 1], (ii) C2, the line 
connecting z = ltoz = l + ioo, (iii) C3, the line connect- 
ing z = 1 + ioo to z = ioo, and (iv) C4, the line connecting 
z = 1 + ioo to z = 0. From the form of fiA, we can easily see 
that J Ca = [because lim z _ i . 00 A(z) — > 00] and J c = — f Cl 
[because A(z) = A(l — z)]. This means that for r« 1 and 
in the regime in which gives the dominant contribu- 

tion to <5 a d(l), the error exhibits an exponentially decreasing 
behavior as 



Ml) 



fT06l 



< 



-wJTr/8 



In this regime, for <J a d(l) < e it is sufficient to have 

TZ 81og(l/e) | 81og(7r/2) 
~ TrJr TrJr 
whereby we can estimate the following run time: 

1 



O -log(l/e) 



(147) 



(148) 



(149) 



This result agrees perfectly with the expected O(N) scaling, 
with a logarithmic (rather than inverse) dependence on the er- 
ror e. Figure Q] depicts the adiabatic error vs time, calculated 
by solving the corresponding Schrodinger equation numeri- 
cally. It illustrates the exponential and polynomial regimes. 

Remark. Here we could not use Eq. ( 1109t because Ha has 
poles zq±. 



2. Constant-norm interpolation 



Here for simplicity, we assume 6(t) = r. From Eqs. ( 11341 ) 
and ( 1135b , the gap $2% is 



A(r) = y/r sec[(l - 2r)(p], 



(150) 



with the minimum value min T A = sjr at r m - m = 1/2, and 
Eqs. (|38) and (|39j yield 

||!T(r)||/J= 1, (151) 
\\H(t)\\/J= ^scc 2 [(l-2r)^], (152) 

the latter with the maximum value max T ||ij(r)j| = <p/^/r 
at T max = 0, 1. Thus, according to the traditional adiabatic 
condition d20"i i, in order to have £ a d(l) < £> we should require 



1 



and in turn, 



= O 



erwr 



(153) 



(154) 



which is larger than the expected 0(\J N /M) Grover-like 
scaling 

On the other hand, truncation of the corresponding 
poly(l/T) expansion j97| ) results in 



^ 2cp 
T > — 

in which we used |ii(l)| = 1 2^2 (0) | = 
[Eqs. CEH and (fT35l) l and ip^n/2. Hence, 

r ra „ = 0(1), 



(155) 
(156) 



which of course is incorrect. 

Now we show that a careful treatment of Zi(t), as in 
Sec.[lV] results in an exponential adiabatic error, and gives the 
correct scaling for the run time. Here, we note that h \ (r) = <p 
[Eq. (TToTll. from which 



Jo 

1 — sin(yj) + cos(t^) tan(y)T / ) 
1 + sin(yj) — cos(ip) tan(<^r') 



arctaiihv/T— r 



JTy/7 



dr', 



(157) 




200 



400 600 
JT 



800 



1000 



FIG. 2: (Color online) <5 a d(l) for the constant-norm interpolation 
[Eqs. 0134b and l |135| ll, obtained by numerically solving the corre- 
sponding Schrodinger equation, for r = 0.01. The dashed lines rep- 
resent the exponential fit e~ JT v^/ 2 [Eq. i(T59}] and the polynomial 
fit 2(p{r)/(JT) [Eq. ( fl55t l. 
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where we used the identity arctanh(x) 
r <C 1 [sin(yj) rj 1], this gives rise to 



iln|4±£|. For 

2 I 1 — x I 



\Mr) 



< 



<P 



[tan(^T') 



J /'[tan(7rr72)]- l ^dT' + 0(^ : ). 



(158) 

Further simplification can be obtained by using the identity 

J 1 [tan(7rr')]"' ia! dT / = sech(7ra/2) (for a > 0) and the in- 
equality scch(y) < 2e~ y (for y > 0); hence, 



Ml) £ Tre-^/a + (V¥). 



(159) 



This is the corresponding exponential behavior for the 
constant-norm Hamiltonian interpolation. It implies that for 
<5ad(l) < e it is sufficient to have 



T> 



or equivalently, 



= O 



21og(l/g) 
log(l/e) 



(160) 



(161) 



which is the expected Grover-like 0(^/N/M) scaling l34l 
l36ll . but with a logarithmic dependence on the error. 

Although Eq. ( 1106b may not be applicable to the case of 
constant-norm interpolation (for Ha does not have any sin- 
gularity), we can apply Eq. (11091 instead. In fact, Eq. ( 1150b 
implies that 



±ioo. 



(162) 



Hence, 



Im 



r z o+ 

/ A(z) dz 
Jo 



17 w — ' (163) 



and in turn [Eq. ( 1109H 

Ml) < R(l) e' JT ^l\ 



(164) 



with some non-exponentially decreasing R(l). This results in 
a scaling similar to Eq. (11611 ). 



3. Constant-gap interpolation 



In this case [A(r) = 1], from Eqs. (|38]>, (|39j, Sffi) , and 
(1138b we obtain 

\\H(t)\\/J= |cos[(l-2r)^]|/^ : , (165) 
||H(r)||/J= ^|sin[(l-2r)^]|/Vf. (166) 

Hence, according to the traditional adiabatic condition ( 1201 ). 
in order to have 5^(1) < e, it is sufficient to have 



9 



Jeyfr' 



(167) 



whereby 



O 



(er) ' 



(168) 



which is quadratically larger than the Grover-like 0(\J N /M) 
scaling. 

On the other hand, noting that |±i(l)| = 1^2 (0)| = 
ip/ x/r(l — r) [Eqs. ( 1137b and ( 11381 )1. the truncation of the 
corresponding poly (1/T) expansion after the first term yields 

2ip 

(169) 



T> 



Je' 



or equivalently 



= O 



(170) 



Clearly, Eqs. (1168b and (1170b are not in agreement. 

It is interesting to note that here one can in fact solve the 
Schrodinger equation exactly. The operator Qt{t) [Eq. (l83b l 
satisfies the following equation: 



O = -iWfl, 

in which 

W = TA^HA- $H A A 

(T37j,(B8},(76} JT f COS[(l — 2r)(p] 



1 +CT; 



(171) 



Pay. 

(172) 



Since W is time-independent, integration of Eq. ( 11711 ) is 
straightforward: 



Mr) 



-i/ T W(l)dr' 



3 in(r^ 2 + (JT)y4) 



^ 2 + ( JTf/A 



I j 

cot (tvV + ( JT) 2 /4) 1 + i(pTo~y - —o-z 



(173) 




100 



F1G. 3: (Color online) <5 a d(l) for the constant-gap interpolation 
[Eqs. {137} and ( Tl38b l. Here, r = 0.001, 0.5 (blue and green, respec- 
tively), and the dashed lines (red) show the 2<p(r)/JT envelopes. 
The points Tk, where 5ad(l)| T = 0, are given by Eq. (T77) , 
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Thus, from Eq. ( l85l ) the adiabatic error is exactly 

I sin y/ip 2 + (JT) 2 /4 



<5ad(l) = </> J 



V<^ 2 + {JTfJl 



(174) 



Figure |3] depicts <5 a d(l) for two different values of r. Note 
that for large evolution times T ^ 2tp/ J, we obtain 



Ml) < 



2ip 
If' 



(175) 



This implies that, in the r< 1 limit, in order for <5 a d ( 1 ) < £ it 
is sufficient to have 



= 



(176) 



which is the Grover-like scaling 0{y/N/M). Notice that in 
this limit the adiabatic error behaves inverse-linearly, <5 ac j ( 1 ) < 
2ip/JT, which is in perfect agreement with Eq. J 169b . In 
addition, we observe that there exist Ts less than the above 
limit in in which the adiabatic error can vanish (hence instan- 
taneous full adiabaticity). According to Eq. ( 1174b . we have 
<5ad(l)| T = where 



JT k = 2y / k 2 ir 2 - ip 2 'w 1 27iVfc 2 - 1/4, 



(177) 



for k € N. The existence of such 7\.s is in agreement with 
Ref. lHH. Figure |3] shows <5 a d(l) for two different values of r. 

Remark. As is evident here the adiabatic error does not 
show any exponential behavior. In fact, neither of the methods 
we discussed in subsection llV B 2l is applicable. 



4. General interpolation 

Here, we discuss the behavior of the exponential 
e - jTim[/ o A] f or tne g enera j interpolation we derived in sub- 
section HV C II Our analysis is based on a formal power series 
expansion of x\ (0(t)) and x 2 (0(t)) in terms of r — recall that 
we are interested in the regime rC 1. We further assume that 
9(t) does not depend explicitly on r. 

Consider the following formal expansions: 

xi{6(r)) = /i(0(t)) + 5l (0(r))r Ql +0(r Ql+£l ),(178) 
s a (0(r)) = / 2 (0(t)) +g 2 (9(T))r^+0(r a ^),(H9) 

in which ot\ and a 2 are some nonnegative numbers (to be de- 
termined later), ex, e 2 > 0, and /x, f 2 ^ 0. We notice that the 
linear and constant-gap interpolations (subsections lIV C 2~l and 
IIV C 4b do not admit expansions as in Eqs. ( 1178b and (1179b . 
Equation (l42l hence yields 

A 2 = (/x - f 2 f + 2(/x - h){9ir ai - g 2 r a *) + 4/x/ 2 r 



ai+Q 2 +ei+e2 



)• 



(180) 



On the other hand, Eq. (1124b yields 

X2i{r) = 1 - 2Vrcot(7r0(r)) + 0(r). 



(182) 



The symmetry x 2 (8{t)) = Xi(l — 0(t)) [Eq. (I127H requires 
that 



9i{0(r)) = 92 (1-9(t)), 



(183) 



which in turn implies oti = a 2 = 1/2. Comparing the terms 
with the same powers of r in Eqs. ( 11811 ) and ( 1182b . we con- 
clude that 



and 



h(9(r))=f 2 (0(r))^f(9(r)), 
gi{0(r)) - g 2 (9(T)) 



f{0(r)) 



= 2cot(7r0(r)). 



(184) 



(185) 



After inserting the above relations back into Eq. d 1 80b and 
using Eq. ( 1185b . we obtain 



A(t) = 2/(0(r))csc(7r0(T))Vf + 0(r). 



(186) 



Now we assume that /(0(r)) ^ everywhere, or if there 
exist points at which / vanishes, their contribution to the inte- 
gral Im[J z ° A] is not substantial. Note that the previous con- 
dition is in fact a condition on the norm of the Hamiltonian — 
because from Eq. ( f38l > 



\H\\/J = 2f + 0(VF). 



(187) 



Therefore, provided that for large times the adiabatic error 
asymptotically behaves as in Eq. dl06b or ( 1109b , we obtain 



Ml) ~ i?(l) e - 2 v^ JT | Im [/o /(e^))c S c(^( z )) dz]\^ 



(188) 



where R(l) is a non-exponential function of T (which may 
also depend weakly on r). Hence, in the r <C 1 limit for 
£ad(l) < £> it is sufficient to have 



T> 



log(l/e) 



2JV? 



Im 



<o(r) 



f(9(z)) CSC(7T0(Z)) di 



(189) 



This in turn implies the following scaling for the run time: 



O 



log(l/e) 



Im 



2 oM 



f(6(z))csc(ir8(z)) dz 



(190) 



„. ., , „ rr=7oi . rn^rn. • / The overall r-dependence here comes from Jr and zn(r); 

Similarly, inserting Eqs. ( 1178b and (1179b into x 2 i = x 2 jx\ r . — ; — v . . w 



gives 



X 2 1 = -p 
h 



h hg\ 



.92 



/f 



^r Q2 + 0( 
Ji 



,ai+a 2 +ei+e2 



)■ (181) 



e.g., we recover the Grover-like 0(y/N/M) scaling if zq does 
not depend on r. This analysis then highlights in a fairly gen- 
eral way the interplay between r, T, and <5 a d(l) m the quantum 
search algorithm. 
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E. A strategy for reducing the adiabatic error 

For most applications it is desirable to make the adiabatic 
error as small as possible. We have seen that 5 a( j(l) contains 
exponential terms, suppressed by the polynomially-decaying 
terms. Therefore, it is useful to somehow extend the domi- 
nance of the exponential term by reducing the contribution of 
the polynomial term, e.g., by prolonging the dominance of the 
exponential regime or by enforcing higher order polynomial 
behavior. In the following we shall discuss a control strategy 
for reducing the adiabatic error by manipulation of the bound- 
ary conditions (see, e.g., Refs. i28[l32"ll ). 



1. A general strategy: Control via boundary conditions 

Equation d97| l demonstrates explicitly how the adiabatic er- 
ror depends on the boundary conditions, up to 0(1/T), Inter- 
estingly, if we choose x 2 (0) = X\ (1) = 0, the prefactor of the 
1/T vanishes, whence <5 a d(l) = 0(1 /T 2 ). In a similar fash- 
ion, one can see that by enforcing suitable (extra) boundary 
conditions on the interpolation x(r) the prefactor of the 1 /T 2 
or even higher order terms can be made zero. This implies 
that by manipulating the boundary conditions of the interpo- 
lation, one may achieve smaller adiabatic errors. This obser- 
vation is a manifestation of the following general theorem: 
If the Hamiltonian H(t) is sufficiently differentiable, forcing 
all time derivatives of the Hamiltonian up to some order k to 
vanish at the boundaries, 



d l H(r) 



r£{0,l} 



= VI e {!,...,&}, 



(191) 



is sufficient for (5 ad (l) = 0(1/T k+1 ) H \M M, \5M ■ It 
is interesting to note that the very same condition together 
with the assumption of the analyticity of H(t) in a small 
strip around the real axis in the complex r-plane give rise to 
Ml) = 0(e- c « JT ), where c(r) = A^ in /max T \\H(t)\\ 2 
[up to an O(l) prefactor] JH. 2 

This is a remarkable result, in that it guarantees that with 
sufficient smoothness at two points one can substantially sup- 
press the adiabatic error. This is a fairly low price to pay for 
higher accuracy. In particular, in experimental realizations, 
manipulating Hamiltonian interpolations only at the begin- 
ning and the end (as opposed to from beginning to end), may 
offer a less demanding control strategy than one seeking to 
control the dynamics instantaneously along the entire evolu- 
tion. 

Now we show that in the framework we developed earlier, 
enforcing the required smoothness properties can be achieved 
by choosing an appropriate 8(t) function. We recall that this 
function was fairly arbitrary; we required that it be a monoton- 
ically increasing differentiable function (8(t) > 0) satisfying 



2 Note that the conditions ±2(0) = 2:1(1) = we obtained above are in 
fact weaker than requiring i'i({0, 1}) = i'2({0, 1}) = [sufficient for 

H ( T )\re{0,l} 



the boundary conditions 6(0) = and 8(1) = 1. We require 
further that 6(t) = $k (t) (for a given k £ N) have the follow- 
ing property: 



d l 8 k (r) 



dr l 



t€{0,1} 



Vie {!,..., k}, 



(192) 



namely, the first k derivatives should vanish at the bound- 
aries. This property is sufficient for fulfilling Eq. d 1 9 1 b be- 
cause H = QdeH. An example of such 8 k (r) is the regular- 
ized incomplete beta function, 



-Xr) 



B T (1 +k, 1 + k) 
Bi(l + fe,l + fc)' 



(193) 



in which B T (a,b) = /JV^l - yf^dy, with 
Re(a),Re(6) > 0, and |r| < 1 1313. 



2. Polynomial reduction 

To demonstrate explicitly how the conditions ( 11911 ) af- 
fect the adiabatic error, we employ the method developed 
in Refs. l28l l32tl for approximating the coefficients in the 
poly (1/T) expansion of 5 a d ( 1 ) . One can construct an 
approximate (unnormalized) ansatz for the solution to the 
Schrodinger equation (O in the powers of 1/T as follows: 



|* n (r)> = e~ iJT f TE ° 



l*o(r))+X) 



1 



1=1 



(JT) 



■Mr)) 



(JT) 



ra+l 



(194) 



with the error 

llllMl)) 



l*n(l)>ll < 



=r||fel(T)>|| 



(JT)n+l 

Here, {\ipi(r))} and {|V'/ L ( T ))} are given as follows [f32h : 

hfo) - a,|*o} + M X 
i 

ai = - J ($oK'r) dr '. «o = 1, 

G r = itH-Eo)- 1 ^ - |*o>(*o|) 



(195) 



Wr) 



(196) 
(197) 

(198) 

(199) 



It is evident that \ip T (r)) = |*oo(r))/||*oo(r)||. 

Provided that all n derivatives of H(t) vanish at the 
boundaries, Eq. ( |191t , all the terms except |$o(l)) an d 
|V>,i + i(l))/( t /T)™+ 1 will vanish JH (see also Ref. jH, 
where with a condition different than Eq. J 198b all the terms 
except |$o) vanish). 

Let us define 



Mi) = Vi - |(Vt(i)|*4i))I 2 /II*4i)II 2 , (200) 



0]. 



S 2 (l) = VI - |(*„(l)|$_(l))| 2 /ll*n(l)|| 2 . (201) 
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Since S(a,b) = yl — |(a|6)| 2 (for normalized |o) and |6)) 
is a distance [Eq. (0], from the triangle inequality S(a, b) < 
S(a, c) + S(b, c) we have: 



M 1 ) <*l(l) + *2(l). 
Notice that <5i (1 ) can also be written as follows: 



(202) 



1(^(1)1^(1))^ 

ll^(l)H 2 ||$oo(l)|| 



1 



o 



|(vl/ n (l)|(|vl/ n (l)) + |^(l)))|^ 
[1 + 0(1/T)] 2 [1 + 0(1/T)] 2 

l-(l + 2Re[(§ n (l)|i? n (l)}]) 
1 



where 



\R n (l))^ |*oo(l)> -!*«(!))• 



(203) 



(204) 



On the other hand, a straightforward calculation (supple- 
mented with induction) shows that 

= 'v/i-K^d)^^!))! 2 /!!^!!)!! 2 



i 2 (o)e +1) (o)| 



yV(l-r)|l-2r| 

+ii 1 (i)e +i) (i)i)+o(^ 



(205) 



This bound holds true for any interpolating paths x(r) for 
which xi(l — r) = X2(t). As can be seen, 6\(1) is negli- 
gible in comparison to <5 2 (1), hence we obtain 



Ml) < 



yV(l-r)|l-2r| 
(JT)™+ 1 

+ii 2 (o)e +i) (o) 



1 



O 



J<n+2 



(206) 



This result is a generalization of Eq. d97| l. For example, in the 
case of the constant-gap interpolation (subsection llV C 4t . this 
error reduces to 



^ad(l) £ 



2ip\ 



,(«+!) 



(1)1 



(JT) 



n+l 



(207) 



in comparison with Eq. ( 1175b . Figure [4] depicts <5 a d(l) f° r the 
constant-norm interpolation [Eqs. (11341 ) and dl35H . It can be 
seen that by increasing k the exponential regime dominates 
longer, while the polynomial regime is pushed farther away to 
the region of large evolution times. However, this improve- 
ment comes at a price. The rate of exponentiality decreases 
with increasing k; that is, if fci < fc 2 and <5 a d(l) < £ for an e 
in the exponential regime for the larger k, then T\(e) > T^e) 
[see the inset of Fig. [4). in other words, for some values of 
e, increasing k might give rise to an increased run time. Of 
course, if e is such that the polynomial regimes dominate for 
both values of k\ and k% the interpolation with the larger k 
(/c 2 ) results in a smaller run time. 



3. Exponential reduction 

Since choosing a^fr) with a larger k benefits the accuracy 
of the adiabatic evolution, it is natural to investigate cases with 
k = oo. An example of such 0oo(t) is 



J Q T b af! (r')dT' 
J'b a0 (r') dr'' 



where 



b a p{r) = e-WV-rW O< a ,0<l, 



(208) 



(209) 



is a symmetric "bump function". We notice that b a p(r) is 
compactly-supported and infinitely differentiable in r € [0, 1] 
(the "Schwartz class" ll68ll ): in particular, 



At 1 



T-e{o,i} 



= G N. 



(210) 



However, it is not an analytic function of r, which implies 
that H (x(0oo(t))) is not analytic either. Despite this infinite 
smoothness, the very lack of analyticity in fact prevents the 
adiabatic error from being identically zero 128112911 . 

Nevertheless, Ooo (r) helps remove the polynomial terms ar- 
bitrarily, hence extending the exponential regime farther. Ad- 
ditionally, the exponent of the exponential term is control- 
lable through varying the parameters a and f3. In this case, 
asymptotic evaluation of the integral d82l with the stationary 
phase method results in a faster-than-polynomial convergence 
to zero. For example, with a = 1 and f3 <C 1 for the constant- 
norm interpolation, we can approach an exponential conver- 
gence similarly to the case with k — 0. The achievement of 
an exponentially small <5 ac j(l) for such functions may be un- 
derstood as an instance of rigorously derived exponentially 
small errors for a general class of functions called the Gevrey 
class iSSHl. 



100F 




FIG. 4: (Color online) 5 a d(l) obtained by numerically solving the 
Schrodinger equation corresponding to the constant-norm interpo- 
lation with Ok(r) the regularized beta function [Eq. d 193t l. for 
k G {0, 1,2} (here r = 2 -8 ). The dashed (red) lines represent 
g-^rJTv^AV) and2v2|0 1 2) (l)|/(JT) 2 [Eq. <207l l. 
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It is evident that manipulating the 0(r) function may also 
affect the exponent lm[]*° A(0(T'))d-r']. This can be ob- 
served in Fig. @] (the inset) through the change of the slope 
of the exponential lines in the log plot. A side consequence 
of using a 0(r) with the desired boundary conditions ( 1192b 
is that increasing k may adversely increase the value of T 
for which <5 a( j ( 1 ) < £ (in the exponential regime). To quan- 
tify how choosing a Ok (r) for the Hamiltonian interpolation 
affects the performance of the algorithm, we propose the fol- 
lowing measure: 



whereby 



<Wi) 



(211) 



namely, the average adiabatic error up to time T, for a given 
k. A larger average error may be interpreted as less efficient 
performance. 

The above problem with the effect of k on exponentiality 
may be partially alleviated in some cases. Recall that in sub- 
section II V C 1 1 we found a fairly general interpolation, which 
resulted in the parametrization x\2(t) [Eq. (I124H — it was 
later that we added further conditions so as to find xi(r) and 
X2 (t) separately. Rather than assuming the condition (1128b , 
let us impose 



2/(0(r)) =1 + C0(t), 



(212) 



for some ( ^ 0. We should be mindful of the fact that, 
from Eq. (1187b . modifying / leads to a modification of the 
norm of the Hamiltonian, or equivalently the maximum en- 
ergy of the system; this is a cost, which should be taken 
care of in the correct estimation of r lun [Eq. d26|il. Addi- 
tionally, from Eq. ( 1186b it is seen that — since we assumed 
/ 7^ — (up to the leading order in r) the gap A vanishes 
where sin(0(zo)) oc oo; i.e., 6{zq±) = ±ioo. In this case, the 
exponent lm[J z " A(0(t')) dr'] [Eq. fT07l l becomes 
r z o + 

(1 + (6(t')) csc(tt0(t')) dr' = 



Im 



Im 



CSC 



>0(t')) dr' 



(213) 



Therefore, by appropriately choosing £ — subject to the condi- 
tion 1 + (8 > — we can tune the exponent of the exponential 
term in 5 a( j(l). This in turn gives us control over the run time 
in the exponential regime. Furthermore, as we argued ear- 
lier, replacing 0(t) — > 6k{r) causes the polynomial terms of 
(5 a d(l) to be 0(1/T k+1 ). Thus, we are now in possession of 
two control parameters k [more precisely 0&(t)] and £ with 
which we can manipulate how the adiabatic error behaves in 
either of the exponential and polynomial regimes. This type 
of control may have applications in experiments in which adi- 
abaticity plays a role. 

As an example, let Ok (r) be the regularized beta function 
[Eq. ( 1193H . For k S> 1 and after employing Stirl ing's approx- 
imation for the factorial function (fc! rj \/27rfc(fc/e) fc 065I1 ). 
we obtain 



niax0(r) = 0(1/2) 



2(2 + fc)^- fe (! + fc)3+fc 



■ Vie, 
(214) 



|F(r)||/J«l + Vfc|a 



(215) 



As a result, for example, in order to keep the maximum energy 
constant, while having the advantages of Ok, we should choose 

lObl = o(i/Vk). 



V. SUMMARY AND CONCLUSIONS 

Adiabatic evolution is characterized by a tradeoff between 
the total time taken and the error in the final state reached, 
relative to the desired adiabatic state. Motivated by a desire 
to understand and optimize this tradeoff, in this work we per- 
formed a detailed analysis of the adiabatic error for the case 
of an adiabatic quantum search algorithm. Rather than using 
the traditional adiabatic condition, with its associated pitfalls, 
we chose to calculate the adiabatic error directly by solving 
the Schrodinger equation. This enabled us to derive an exact 
relation for the adiabatic error. Building on this exact result, 
we employed a formal polynomial series expansion in 1/T 
for calculating the error term by term. This also allowed us 
to bound the adiabatic error. We showed that the polynomial 
expansion should be truncated carefully if we aim to obtain a 
reliable estimate for the run time of the algorithm. 

We demonstrated that employing a different technique 
based on complex analysis reveals, in fairly general situa- 
tions, a regime of exponential time-dependent decay of the 
adiabatic error, preceding a polynomial regime. The latter 
has been shown to be a general feature of adiabatic Hamil- 
tonians, whereas the existence of an exponential precursor is 
not always guaranteed. We showed how, in case these two 
regimes both exist, one can estimate the crossover region — 
i.e., the time at which a transition between the two regimes 
takes place. Equipped with this, we provided an estimate 
for the minimum time required for the algorithm to achieve a 
given accuracy threshold. Discerning the exponential regime 
enabled us to give an improved total time estimate, circum- 
venting the overestimate arising from the error bound using 
only the polynomial expansion. Indeed, the total time esti- 
mated from the exponential regime always gave the correct 
scaling with system size (the well-known quadratic speedup 
over classical search), while the estimate resulting from the 
polynomial regime resulted in unreliable and even erroneous 
results. 

We also obtained a specific class of Hamiltonian interpo- 
lations for the search problem. To this end, we employed 
a recently developed theory, based on the geometry of adi- 
abatic evolutions, for obtaining suitable adiabatic interpola- 
tions. This theory separates the adiabatic error into geometric 
and non-geometric parts, and minimizes the former. We dis- 
cussed three special cases of the resulting class of interpola- 
tions in detail: (i) linear, (ii) constant-norm, and (iii) constant- 
gap interpolations. 

Finding strategies for minimizing the required total time 
as a function of a given resource (system size, for exam- 
ple) is a desirable goal for many applications, and is also of 
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fundamental importance for the control of quantum systems. 
We demonstrated explicitly how by imposing fairly general 
controllability assumptions, which should be experimentally 
straightforward to realize in certain scenarios, one can achieve 
a significant reduction of the adiabatic error, and hence im- 
prove algorithmic performance. The method we used relied 
on a polynomial expansion of the adiabatic error, and resulted 
in the suppression of polynomial terms in 1/T by requiring 
smoothness for interpolations at the initial and final times. It 
is evident that controlling the interpolation in this manner, at 
only two points, has a substantial advantage over instanta- 
neous control of the Hamiltonian along the entire evolution. 
However, we demonstrated that there is an extra price to pay 
for this error reduction: the exponential regime (if exists) is 
extended, but with a slower rate of decay. This, in turn, may 
result in an overestimation of the run time of the algorithm for 
some values of the error threshold. We proposed a measure for 
quantifying the performance of an adiabatic interpolation with 
various controllability properties. In some cases, we also sug- 



gested a remedy for the above problem. This fix necessitated 
further control over the Hamiltonian interpolation, directly re- 
lated to the amount of accessible energy in the system. The 
interplay between the degree of required control over Hamil- 
tonian interpolations and the run time needed for achieving a 
given accuracy was thus clearly exhibited. 

Although we focused on the quantum search problem, our 
methods and most of our results are applicable (perhaps with 
minor modifications) to a wider class of problems — as dis- 
cussed in the text. Since a principal goal in adiabatic quantum 
algorithms, adiabatic quantum transport, quantum annealing, 
and other applications of the adiabatic theorem, is the design 
of algorithms with favorable performance-resource tradeoff, 
we hope that our results will be of use in related physical ap- 
plications. 
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